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Abstract 

The squashing factor Q (Titov et al. 2002), a property of the magnetic field line mapping, 
has been suggested as an indicator for the formation of current sheets, and subsequently 
magnetic reconnection, in astrophysical plasmas. Here we test this hypothesis for a particular 
class of braided magnetic fields which serve as a model for solar coronal loops. We explore 
the relationship between Quasi-Separatrix Layers (QSLs), that is, layer-like structures with 
high Q value, electric currents and integrated parallel currents, the latter being a quantity 
closely related to the reconnection rate. It is found that as the degree of braiding of the 
magnetic field is increased the maximum values of Q increase exponentially. At the same 
time, the distribution of Q becomes increasingly filamentary, with the width of the high-Q 
layers exponentially decreasing. This is accompanied by an increase in the number of layers 
so that as the field is increasingly braided the volume becomes occupied by a myriad of thin 
QSLs. QSLs are not found to be good predictors of current features in this class of braided 
fields. Indeed, despite the presence of multiple QSLs, the current associated with the field 
remains smooth and large-scale under ideal relaxation; the field dynamically adjusts to a 
smooth equilibrium. Regions of high Q are found to be better related to regions of high 
integrated parallel current than to actual current sheets. 
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1 Introduction 



Several different mechanisms have been proposed as explanations for solar coronal heating and 
indeed it is likely that a number of different heating processes are responsible for the extreme tem- 
perature values found. Much of the corona is considered to be close to force-free, that is, in a state 
satisfying j x B ~ 0. This consideration led to the proposal of one of the earliest and most debated 
coronal heating theories, that of topological dissipation (Parker, 1972). Parker theorized that the 
coronal magnetic field is unable to relax to a smooth force-free equilibrium following arbitrary per- 
turbations of the field via footpoint motions and the consequence is a development of tangential 
discontinuities in the field. Magnetic rcconnection may then occur across the resulting current 
sheets. Since Parker's original paper a multitude of arguments for and against the hypothesis have 
been given but no unanimity presently exists (e.g. van Ballegooijcn, 1985; Galsgaard & Nordlund 
1996, Craid & Sneyd 2005). 

One technique presently employed for predicting current sheet formation in continuous fields 
is to examine a quantity known as the squashing factor, Q (Titov et al. 2002; see Appendix 
A for further details of the function itself). This is a measure for continuous magnetic fields 
which describes the level of squashing of an infinitesimal flux tube; Q is large in regions where 
the footpoint mapping of the field is highly distorted. Regions of very high Q outline so-called 
quasi-separatrix layers or QSLs (Priest & Dcmoulin 1995) and can be considered as an analogue 
of the separatrix surface in configurations with no null points. QSLs are thought to be sites of 
preferential current sheet formation (Demoulin et al. 1996; Titov et al. 2002) and a number of 
numerical experiments have supported this hypothesis (e.g. Longcope & Strauss 1994; Milano et 
al. 1999; Galsgaard et al. 2003; Aulanier et al. 2005). In addition, observational studies have 
attempted to examine the relation between QSLs and current sheets (e.g. Schmieder et al. 1997; 
Fletcher et al. 2001) although this is difficult due to the spatial resolution and quality of the data 
required for an accurate depiction of the QSLs. 

In an earlier paper (Wilmot-Smith et al. 2009) we considered a class of braided magnetic fields. 
The braids were initially given in closed form (and not in MHD equilibrium), then inserted into 
a numerical scheme and underwent magnetic relaxation towards a force-free equilibrium. In that 
paper the current structure of the fields was looked at in some detail. The results are summarized 
for convenience and put into context in Section 3. Here we return to these models for braided fields 
in order to examine the relationship between QSLs, currents and a quantity which is important 
for magnetic reconnection in three dimensions, the electric current integrated along magnetic field 
lines. As such we begin by describing the braided fields under consideration. 

2 Model for a class of braided fields 

In Wilmot-Smith et al. 2009, (hereafter Paper I) we introduced a class of braided fields designed 
to model magnetic fields in solar coronal loops. The fields are given as E n , n e Z + , where E n 
can be considered as a concatenation of n times the 'elementary' field E. The braided field E 3 is 
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modelled on the pigtail braid as illustrated in Figure [T] 




E 



(a) (b) 



Figure 1: (a) The pigtail braid and (b) some field lines from the field E 3 modelled on the pigtail 
braid. E 3 may be considered a concatenation of three elementary parts, E (indicated in the figure) 
and, accordingly, E n may be constructed for various n via similar concatenations. 

The nature of E and the construction mechanism for E n may be seen from that figure. The 
field E consists of two localized regions of twist in an otherwise uniform field. These two regions 
are of opposite sign but are of the same underlying form - the superposition onto a uniform vertical 
field, boi, of a toroidal flux ring given by 

B e = 26d*£«p 

In Cartesian coordinates, placing the centre of this region of toroidal flux at the point (x c ,y c , z c ) 
we have 

T3 o b ° k ( -(x-x c ) 2 -{y-y a ) 2 (z-z c ) 2 \ 

B c = 2 exp — (- (y - y c ) x + (x - x c ) y , 1 

a \ a 2 1 ) 
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where c = {x c , y c , z c , k, a, 1} and bo = 1 is taken throughout. The closed form expression for E is 
given by 

2 

b z + ]T B Ci 

»=i 

where C\ = {1,0, -4, 1, y/2, 2} and c 2 = {— 1, 0, 4, -1, y/2, 2} and, similarly, that for E n by 

2n 
i=l 

with appropriate choices for the parameters Cj. For each field £ n we are able to find the equations 
of the magnetic field lines in closed form and hence an exact expression for Q for each n can be 
easily calculated. 

3 Current structures 

In paper I a Lagrangian numerical scheme (Craig & Sneyd, 1986) was used to carry out an ideal 
relaxation of E 3 towards a force-free state with the field on the boundaries of the domain fixed. 
The maximum Lorcntz force within the domain decreased by two orders of magnitude during the 
relaxation and the resultant field is close to force-free (for a precise definition of what is meant by 
'close' see Section 3.2 of Paper I). Beyond this point problems with the numerical accuracy of the 
scheme (Pontin et al, 2009) mean no further relaxation is possible until the numerical difficulties 
are resolved. In the relaxed state for E 3 the current is smooth and of a large-scale (see Figure 3 
of Paper I) - two regions of current extend vertically throughout the domain and the maximum 
magnitude of the current, (j m ax = 1-47), is around half of that of the initial state. No evidence for 
current sheet formation in the ideal relaxation was found. 

It is well known that continuous motions on the boundary of an initially smooth magnetic field 
cannot lead to truely singular current sheets under an ideal (line-tied) evolution (Van Ballegooijcn 
1985). However it might seem that the lack of current sheets in our relaxed state contradicts 
previous results, for example those of Longbottom et al. (1998) and Galsgaard et al. (2003). These 
investigations found the formation current concentrations under shearing motions on the boundary 
of an initially current-free configuration. Since we can think of our braid as having been created 
by a sequence of such shearing motions applied to an initially homogenous field, we might also 
expect to find similar current concentrations in the relaxed state. A closer comparison, however, 
shows that the twist, measured by the ratio of the displacement distance (Id) to the length of the 
domain (L) in which current concentration are observed is comparatively high in the experiments of 
Longbottom et al. (1998) and Galsgaard et al (2003). In particular, Longbottom et al. (1998) looked 
in detail into this question and found significant build-up of current for values above ld/L w 0.6 
(ld= 'shear amplitude' in their terminology and L = 1). If we were to create our elementary braid E 
(which comes closest to the configuration they considered) by a similar shearing pattern, we would 
obtain a ratio ld/L < 0.18 (Id ~ 2\/2; L — 16). This is a value for which no noticeable current 
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concentration has been observed (see Figure 2 of Longbottom et al. 1998) in agreement with our 
findings. A similar argument applies to the conditions under which Galsgaard et al. (2003) found 
current concentrations, although their experiment does not allow a direct comparison as they used 
a resistive code. In their case, values of Id ( 'shear distance' in their terminology, L = 1) below 0.2 
also showed no significant build up of current. 

While this shows that our results do not contradict previous findings of current sheet forma- 
tion, it is perhaps still surprising that a smooth near force-free equilibrium for such a complex 
configuration exists. In order to understand this one has to bear in mind that our equilibrium 
is not an exact force-free state but only close to a force-free state (for an exact definition how 
close see Section 3.2 of paper I). If the residual forces in the relaxed state could be balanced by 
a pressure, our relaxed state would have a plasma beta of (3 « 0.009; although this is well within 
the quality of 'forcc-frccness' usually assumed for magnetic loops in the solar corona, it is not zero. 
Moreover, the field has comparatively low maximum force-free field parameter a (when measured 
in a dimensionless way) and it is known that for sufficiently small values of this quantity, force- free 
equilibria exist for arbitrarily braided fields (Bineau, 1972). Unfortunately the proof given does 
not give an explicit upper bound but only provides the existence of such a number. These two 
properties, a low a and an only approximately force-free state, appear to allow for a sufficiently 
large space of smooth solutions to encompass the braids we have investigated. 

Whilst in two-dimensions the current itself is the crucial quantity for magnetic reconnection, in 
three-dimensions the integrated parallel electric field along magnetic field lines plays an important 
role (Hesse & Schindler, 1988). In resistive MHD, the parallel electric field is related to the 
parallel current by the relation J E\\dl = r] J j\\dl (assuming a uniform resistivity -q). Motivated 
by this consideration, in Paper I we examined the integrated parallel current, hereafter denoted 
by J7|| , in both the initial and relaxed states for E 3 . In both cases a highly filamentary structure 
containing multiple thin J|| layers was found. During the relaxation the spatial structure of J\\ 
was approximately conserved. More generally, by considering E n for various n, the thickness, d, 
of the J|| layers was found to decrease exponentially with n (degree of braiding) . 

A quantity proposed as an indicator for the formation of current sheets (and hence for the 
occurrence of magnetic reconnection) in continuous fields is known as the squashing factor, Q 
(Titov et al. 2002; Titov et al. 2009). Layers in which Q is large, so-called quasi-separatrix layers 
(QSL), can be understood as a generalisation of separatrix surfaces for cases where no null points 
exist in the domain under consideration. Our previous findings motivate a consideration of the 
nature of Q in the class of braided magnetic fields E n . We relate these findings both to the current 
structure and the integrated parallel current structure of the fields. 

4 Squashing factor and QSLs 

A method for obtaining the squashing factor Q in the braided fields considered is outlined Appendix 
A. Note that for E n the squashing factor Q and the newer definition of the slip-squashing factor 
(Titov 2007, Titov et al. 2009) give the same result since the magnetic field is normal to both 
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Figure 2: (left) Q as a function of field line positions on the lower boundary of E and (right) 
integrated parallel current along field lines, on the same boundary for the field E. 



z-boundaries. 

We begin the discussion by considering the simplest field, E. Under ideal relaxation towards 
a force- free equilibrium, the current, j, in E remains smooth and large-scale, with the maximum 
value being j max — 1.07 in the relaxed state (compared with \ max = 2.83 initially). The squashing 
factor Q for E is shown in the left-hand image of Figure [2] as a function of field line positions on the 
lower boundary. High values of Q (regions where the field line mapping is strongly distorted) are 
found in two layer- like regions, and the maximum value of Q is Q\, ma x = 241.5. Contours of the 
parallel current integrated along field lines, J7j|, are shown in the right-hand image of Figure [2j The 
characteristic scales of J\\ are larger than those of Q and no relation between the two is evident. 




Figure 3: (left) log 10 (Q) as a function of field line position on the lower boundary and (right) 
integrated parallel current along field lines, \J\\\, on the same boundary for the braided field E 3 
(based on the pigtail braid) . 
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Next consider the field E 3 , that based on the pigtail braid and on which most of the analysis 
of Paper I was based. Recall that a smooth equilibrium with only large-scale current structures 
was obtained in an ideal relaxation of this field. 

The squashing factor Q on the lower boundary of the domain is shown in the left-hand image 
of Figure [3] Here the maximum value of Q is Q3 <max = 7.17 x 10 5 (with Q calculated at 1396 2 
points on the lower boundary) and regions of enhanced Q occur in a multitude of thin layers in 
a large portion of the domain. Between these layers Q drops significantly, even to its minimum 
possible value of 2 in many locations. A typical half-width at half-maximum (HWHM) of the QSLs 
is ~ 10~ 3 or 0.01% of the domain width. 

The integrated parallel current structure, J'n, for E 3 is shown on the lower boundary of E 3 in 
the right-hand image of Figure |3] A resemblance between Q and J\\ is shown, with both quantities 
having a similar global structure and containing thin layers of approximately equal width. However, 
Q is more filamentary than J7|| and has a much greater range of values. In order to make a more 
precise comparison, the reader is referred to the upper panel of Figure [3] where both quantities are 
shown for only a subsection of the domain. 

Finally consider the field E e . Both Q and J\\ are shown over a subsection of the domain in 
the lower panel of Figure [4] The structures of the two quantities are highly filamentary, both 
containing extremely small scales (HWHM of ~ 10~ 5 ). It is seen that Q and J\\ here have a 
strong resemblance in terms of the global pattern of their thin layer-like structures. However, the 
locations of the regions of high (low) Q and high (low) Jn do not exactly coincide. 

Considering now Q alone for the various fields E n , we seek to determine how the maximum 
value of Q and the characteristic width (specifically HWHM) , d,Q , of the layers of high Q depend 
on n. It is found that the maximum of Q for E n , Q n ,max, in the domain increases exponentially 
with n - as shown in Figure^ a) - for which a best fit line is 

Qn,max — 2.97 X 10 

In Paper I a linear increase in the maximum of J\\ with n was found, as expected since the vertical 
length of the domain is also increasing linearly by construction of the fields. The typical width, 
<1q, of the regions of enhanced Q decreases with n. More precisely, there is an exponential decrease 
in cLq with n - as demonstrated in Figure [5jb) - for which a best fit line is 

d Q = 1.90 x 10~ - 93 ". 

In Paper I a similar exponential decrease in the width (dj^ ) of regions of high J7j| was found. 

We now address the dependence of the maximum value of Q, Q n ,max on n. Recall that Q n ,m,ax 
was found to increase exponentially with n. Notice that by the construction of E n we may place 
an upper bound on the maximum value of Q from the method of its calculation by composition of 
matrices. Letting Qi denote the maximum value of Q for E l , then the submultiplicative property 
of the p = 2 norm tells us that Q n ,max < (Qi,max) n ■ Evaluating, we obtain the upper bound 
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Figure 4: Figure showing log 10 (Q) and the absolute integrated parallel current along field lines J\\ 
for E 3 (top panel) and E 6 (lower panel). For clarity, only a subsection of the domain has been 
shown. 
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Figure 5: (a) Maximum value of Q with n in the full domain. An exponential exponential increase 
in Q with n is shown, (b) HWHM of Q layers with n. An exponential decrease in Q with n is 
found. Here the width has been taken as the shortest such width along a cross section of the 
domain on the lower boundary. 
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as (Qi,max) n ~ 10 4 7n /200, whilst the true dependence is given by Q n ,max ~ 2.97 x 10 18n . The 
dependence on the smaller power then comes from the filling factor of the field E, i.e. the number 
of points on the lower boundary of E with highest Q\ that get mapped under F to points with 
high Qi for the next composition (for E 2 etc.) 

In the next section we seek to explain these findings and explore theoretically the relationship 
between Q and the current structures. 

5 Relating QSLs and current structures. 

As mentioned above, the squashing factor has been suggested as an indicator for the formation 
of current sheets. Layers in which Q is large outline quasi-separatrix layers (QSLs). The link 
between QSLs and current sheets relies on the following hypothesis (referred to below as the 
QSL-hypothesis). Where Q is high there is a sensitive dependence of the end points of magnetic 
field lines on the starting points (or vice- versa since Q is symmetric). Hence certain motions 
of the plasma on one boundary (particularly those crossing the QSLs) would lead to very high 
velocities at the other boundary, provided the evolution is ideal and the magnetic configuration is 
unaffected by this perturbation. For sufficiently high Q this is inconsistent with v < and so 
the magnetic configuration has to change. Then a lack of neighbouring smooth equilibria can lead 
to the formation of current concentrations or even current sheets. 

Within the experiments described above we are in a position to test the QSL-hypothesis, 
given that we can think of the relaxed configuration for E 3 as as an ideal deformation of the 
initial configuration. The initial configuration shows an abundance QSLs but no strong current 
concentrations. In the subsequent ideal relaxation the QSLs remain but no current sheets (in the 
sense of a sheet-like structure of j) form. This result appears to contradict the QSL-hypothesis. 
There arc, however, two ways in which the system can escape the necessity of forming current 
concentrations . 

Firstly, the deformation applied could be unfavourable for the formation of current sheets 
(Titov et al. 2003). For example, Galsgaard et al. (2003) considered a configuration containing 
two intersection QSLs and applied two types of deformations on the boundary. They found that 
only one of the types of deformation lead to significant build-up of current. 

Given the complexity of our configuration, it seems unlikely that all of the motions which occur 
during the relaxation are of the type unfavorable for the formation of current sheets unless the 
nature of the deformation applied (v~ jxB) excludes any such favorable deformation. The second 
possibility is that following the applied deformations the system is able to adjust to a neighbouring 
smooth equilibrium. This possibility is supported by the arguments given in Section 3 that our 
final state is a low alpha, near force-free state. 

Although current layers are not found in our experiment, it does show layer-like structures of 
. Such layers have a broader definition in that only the integral of j and not j itself needs to 
show small scales. The following example, however, shows that there is no simple relation between 
cither Q and j or between Q and J\\ . We have chosen this example because it is locally similar to 
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Figure 6: Q (solid line) and |J7j|| (dashed line) for the illustrative example of a particular ID force- 
free field (described in Section 5). Here \ \ = 4|j|. In this example Q is not a good predictor of 
current features. (Parameters: c\ — 0, c% = 1, A = 1, s = 4.) 

the braided fields discussed, in that in both cases Q and J\\ each have a layer-like structure, with 
multiple layers in which the quantities are enhanced. Consider the force-free magnetic field 

B = (ci sin (sin (Ax) /A) + c 2 cos (sin (Ax) /A)) y + (ci cos (sin (Ax) /A) — c 2 sin (sin (Xx) /A)) z 

for which j = aB with a — cos Xx, and take z € [0, s], x, y, € (— oo, oo). On the lower boundary, 
z = 0, we have 

Q = 2 + s 2 cos 2 (Ax) (ci cos (sin (Ax) /A) — c 2 sin (sin (Ax) /A)) 2 , 
= / j • B <is' = s(c 2 + c 2 ,) cos Ax, 

Js'=0 

and 

|j| 2 -( c 2 + c 2 )cos 2 (Ax). 

Here the relationship between Q and currents crucially depends on the choice of parameters. 
Setting C2 = 0, for example, maxima and minima of Q, |J7||| and |j| coincide. Setting instead 
Cx = 0, then minima in Q correspond to both maximal and minimal regions of |j7j|| or |j| (see 
Figure [6|. Here we see that the maxima of Q increase quadratically with the vertical extent (s) of 
the domain under consideration, whilst the maxima of J7j| increase only linearly and the maxima 
of |j| are independent of s. 

Note that this is a specific illustration of a more general case of a force-free field B — B y (x)y + 
B z (x)i with force-free parameter a(x). Taking z £ [0, s] then on the lower boundary Q = 2 + 
s 2 a 2 B 2 z while j - a (S 2 + B 2 ) 1 ' 2 and / j^dl = sol (El + B 2 ) 1 ' 2 . 

The above described situation corresponds well to the braided fields E n in that both Q and J7j| 
have a layer-like structure. For E n we also find that there is no exact correspondence between the 
layers of high Q and greatest J\\ . However, there is a similarity between Q and J\\ in terms of then- 
global structure, particularly for moderate and large n (see Figures [3] and |4| . This similarity comes 
not from an inherent relation between Q and J7|| but rather from the fact that both quantities are 
functions of the field lines. 
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Figure 7: Illustrative example demonstrating properties of the field line mapping F n (x,y). For 
the function -0 = exp(— x 2 /4 — j/ 2 /4), the left hand image shows ip(F(x,y)) centre image shows 
4>(F 2 (x, y)) and the right hand image shows ip(F 3 (x, y)). As F is repeatedly applied to the function 
the image develops small scales; for this reason since Q and J7j| take (by definition) one value on 
each field line, they too must develop small scales on the boundaries of E n for increasing n. 

To see this, we prescribe a smooth function (ip(xo,yo)) of field lines on the lower boundary of 
each of the braided fields E n , n = 1,2,.... We then map this function (i/j) along field lines of 
E n to the upper boundary, i.e. find ip(F n (xo, yo))- This is demonstrated for a specific example 
in Figure [7] for E l , i — 1, . . . , 3. As F is applied repeatedly to any function ip(xo, yo), then the 
characteristics of the map will result in a stretching and contracting of tp in some regions of the 
domain according to the nature of F at those locations, specifically the Lyapunov exponents of the 
map. After repeated application of F any two functions will appear similar; the apparent structure 
of the two functions will be determined by the mapping F. 

From this observation we deduce that the structure of both Q and J7j| will depend crucially 
on the plane in which they are viewed and, by the symmetry of the braid, have the same scale 
on both the lower and upper boundaries. Accordingly, J7j| will, for sufficiently high n (about 3 
in this case) , have a filamentary distribution on the lower boundary that is reflective of the field 
line mapping rather than the 'underlying' nature of the parallel currents. These properties are 
confirmed by illustrating Q and J\\ for E 3 in a cross-section in the middle of the domain (z = 0), 
thereby subtracting as much as possible the effect of the field line mapping. This is shown in 
Figure [8] where the quantities are seen to attain larger scales than on the domain boundaries and 
less similarity between Q and J7j| is apparent. 

6 Conclusions 

In this paper we have examined the link between the squashing factor Q and current structures in 
a class of braided magnetic fields. The squashing factor (Titov et al. 2002) is a geometric measure 
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Figure 8: (left) log 10 Q as a function of field line position in the central plane (z = 0) for the field 
E 3 and (right) integrated parallel current along field lines, |j7j||, for the same field, also in the 
central plane. 



for continuous fields that serves to identify regions where the mapping of magnetic field lines is 
highly distorted. It is used to identify quasi-separatrix layers (QSLs, Priest & Demoulin 1995) - 
features of continuous fields which provide an analogue to separatrix surfaces in fields containing 
null points. Previously, regions of high Q have been considered to be a reliable indicator for the 
formation of current sheets under certain deformations on the boundaries. 

Here we are interested in both the relationship between Q and the current j, as well as that 
between Q and the parallel current integrated along field lines, J\\. The reason for this is that 
in 3D a key quantity for magnetic reconnection is the parallel electric field integrated along field 
lines. In resistive MHD with uniform resistivity the parallel electric field and parallel current are 
related via 

J E\\ dl = r] J j|| dl = r)J\\ . 

We have considered a particular class of braided magnetic fields with no net twist, the fields 
being labelled E n (n = 1, 2, . . .), where an increase in n results in a field of increased complexity. 
In an earlier paper (Wilmot-Smith et al. 2009, Paper I) the field E 3 was taken as an initial 
condition in a magnetic relaxation simulation. The field on the boundaries was held fixed and a 
magnetofrictional code used to carry out an ideal relaxation toward a force-free equilibrium. In 
that process the spatial scales associated with the current, j, were found to remain large, suggesting 
a smooth equilibrium corresponding to each E n can be found. 

In Paper I the structure of J7j| for E n was examined. This quantity was shown to display small 
spatial scales and to have an increasingly filamentary structure with increasing n. This structure 
was preserved in the ideal relaxation process. The consequence was shown to be that for a coronal 
field with any finite resisitivity, after a certain degree of braiding via photospheric motions the 
field will undergo a loss of equilibrium due to the high gradients in J7j| , regardless of the current 
structure itself. 

Here we have investigated the nature of the squashing factor Q and associated QSLs for the 
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fields E n . A method was given (Appendix A) to calculate Q exactly, a property that becomes 
important in the analysis where extremely small scales in the quantity are found. It was shown 
that for the field with the least structure, E, Q is enhanced in two layers of width ~ 5% of the 
domain width, i.e. there are two quasi-separatrix layers (QSLs) present. As the degree of braiding 
is increased, Q develops an increasingly filamentary structure, with many narrow layers of enhanced 
Q present. These are true layers in the sense that Q drops to low values between them. The typical 
width of the layers decreases exponentially with n - for E 3 a typical layer width is 0.01% of the 
domain width. Additionally, with increasing n the maximum value of Q exponentially increases. 
The result of the braiding process is shown to be that the domain becomes filled with a multitude 
of QSLs. 

Despite these findings, we find that under ideal relaxation the braided fields E n do not develop 
any small-scale structures in the current j. This demonstrates an alternative to the view that, in 
a configuration with QSLs (outlined by regions of high Q) plasma motion is highly likely to lead 
to the formation of current sheets along the QSLs - another possibility is that the magnetic field 
may adjust to a neighbouring smooth equilibrium. This equilibrium may itself still contain QSLs. 

We also considered the relationship between Q and the integrated parallel current, J7j|. In the 
simplest braid, E, the two quantities are unrelated. With increasing braid complexity (n), Q and 
J7ji have some similar characteristics; for moderate and large n both quantities arc filamentary 
with layer-like regions in which high values are obtained. However, no simple relationship has been 
found between the locations of maxima in the two quantities. 
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Appendix A 



Here we give a method for calculating the squashing factor in the braided fields E n . For this 
we must obtain the mapping of field lines from the lower boundary of the domain to the upper 
boundary. We make use of a property of the fields E n , specifically that E n is a concatenation of 
n times the 'basic' field E. 

Consider first the field B c + b^z (see Equation [l]) , that is, a single toroidal flux ring imposed on 
the background uniform field. The equations X (xq, s) of the field line passing through the point 
x o = (xo, 2/o j z o) are obtained by integrating 

where the parameter s is related to the distance A along field lines by dA = \B\ ds. The components 
of X (xo, s) are given by 

X = (x - x c ) cos C - (2/0 - y c ) sin C + x c, 

Y = (yo ~ Vc) cos ( + (x - x c ) sin £ + y c , (2) 
Z = b Q s + z , 



where 

c = fcv^exp ( ^^-^jiyy-y^ ( erf ^ 

For the field E, the two toroidal flux rings have been placed sufficiently far apart compared 
with their characteristic length scales that the field in a plane equidistant b between them can be 
approximated as vertical, boz. Accordingly, we may replace the error functions in Equations Q 
by ±1 and obtain an expression for the mapping of field lines from the point (x, y) in a 2D plane 
(perpendicular to the field) below the single region of twist in B c + b z to the point f (x,y) in 
similar a plane above it: 

f p : R 2 -> W 2 

f p : (x, y) -> ((x - x c ) cos y c ) sin ^ + x c , (y - y c ) cos $ + (x - x c ) sin £ + y c ) , 

where p = {x c , y c , k, a, 1} and 

^ = 2V7rfc a exp[ I. (3) 

The mapping of field lines from the lower to upper boundary of the field E may then be found 
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by composition of f p with itself using with the relevant parameters for each of the mappings: 

F(x,y) = f p2 o f pi! 

where pi — {1,0, 1, 2} and p 2 — {1, 0, —1, \/2, 2}. A similar mapping, F n (x, y), of the field lines 
from the lower to the upper boundary of E n is obtained by repeated composition of F. For example, 
the field line mapping from the lower to upper boundary of E 3 is F 3 (x, y) = F o F 2 = F o F o F. 
The squashing factor Q for E n can be found directly from the Jacobian (DF n ) of the mapping 

pn. 

||DF"|| 2 

Q 



det (DF n ) 

where || ■ || denotes the entry wise p = 2 norm. In other words, letting 



DF" 

say, then 



/ dX 


ax \ 




( a b \ 


1 Ox 


dy 


H 


\ dY 


dY 




V c d) 


\ dx 


dy / 





a 2 + b 2 + c 2 + d 2 



ad — be 

For the fields E n the determinant of the Jacobian mapping DF™ is exactly unity (det(DF™) = ad— 
be = 1) since on both boundaries the field is the constant field b^i. This method determines Q in the 
initial states for E n . Since in the ideal relaxation towards a force-free equilibrium the boundaries 
are held fixed, the structure of Q on the boundaries of the domain is preserved throughout the 
experiments. 

Using the same approximation as above, we are able to find a closed form expression for the 
integrated parallel current along field lines in the initial states for E n . Taking only a single region 
of toroidal twist, the field line with initial location (x,y) in a plane below the region has a total 
integrated parallel current of 



Jn one = f j\\dl = — (a 2 - - x c f -(y- y c ) 2 ) £, 

"'twist / one jj,a V / 



twist 



where £ is defined in Eqn. The total integrated parallel current along field lines for E n may be 
found using repeated composition of this expression in a manner similar to that used for Q. Note 
that in Paper I J7ii was evaluated directly via numerical integration. 
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